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Abstract. - This Letter shows that for particularly shaped background particle distributions 
momentum exchange between phase space holes and the distribution causes acceleration of the 
holes along the magnetic field. In the particular case of a non-symmetric ring distribution (ring 
with loss cone) this acceleration is nonuniform in phase space being weaker at larger perpendicular 
velocities thus causing deformation of the hole in phase space. 
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c/} Introduction. — In configuration space, phase space 
Wholes appear as localised intense electrostatic fields E\\ h = 
^O— Vy^/i with broadband spectral signature parallel to the 
c/^ ambient magnetic field. In velocity space they form nar- 
£^row regions of lacking particles of one signature kept alive 
^^for a limited time by the electrostatic field. Ion holes are 
i_local deficiencies of ions while electron holes are local defi- 
ciencies of electrons. Thus the former correspond to weak 
t 7T negative, the latter to positive space charges Qi, e - In this 
Letter we deal with electron holes which can be excited by 
beam or current instabilities parallel to the ambient mag- 
^^netic field B, like the two- stream instability which works 
'^for electron drifts Vd > v e , larger than the electron ther- 
^^mal velocity v e [3,4]. At lower drifts this instability is 
replaced by a modified version [14] which is a form of the 
modified two-stream instability [13,16,17,24]. Their the- 
T".ory has been given by Bernstein, Greene and Kruskal [2], 
>Schamel [30-32], Dupree [8,9] and Turikov [38]. Simula- 
k^jtions by Newman et al. [22,23], Muschietti et al. [20,21] 
^_jand others have shown that electron holes are the natu- 
C^ral nonlinear state of these instabilities, being Debye scale 
entities along B in configuration space, and of short exten- 
sion in the parallel velocity component v\\ in velocity space. 
They contain a dilute component of trapped electrons of 
density N t of low energy mv\ /2 < | ecj>h | . In configuration 
space they are oblate in the direction perpendicular to the 
magnetic field (pancakes). Their behaviour in v± has not 
yet been investigated in detail. It is, however, reasonable 
to assume that the holes are either gyro-limited, being 



of transverse spatial extension up to the thermal gyrora- 



dius A 



h± ~ r c 



v e /co c or inertia limited A 



h± 
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Their life time is determined by the stability of the holes 
with respect to the generation of whistlers, trapped par- 
ticle instabilities, particle trapping, heating and diffusion 
and the corresponding generation of dissipation (see, e.g., 
Newman et al. [23]). One might believe that these micro- 
scopic entities are of minor importance for the behaviour 
of the plasma. However, in collisionless plasmas they form 
an important dynamical source of dissipation. They heat 
and accelerate electrons, cause beam cooling, and are sus- 
pected to provide a substantial part of the dissipation that 
is needed in collisionless shocks and in reconnection. In 
collisionless shocks they might contribute to the emission 
of radiation causing the badly understood type II bursts. 
Some time ago we proposed [25] that phase space holes 
contribute to electron cyclotron maser emission [35] gen- 
erating auroral kilometric radiation in the upward current 
source region where the holes have been identified [26] 
subsequently, forming what we called 'elementary radi- 
ation sources'. For this to work the holes must become 
deformed in phase space in order to attain a perpendicu- 
lar phase space gradient dF(y\\ , v±)/dv± on the electron 
distribution function, which is required by the cyclotron 
maser mechanism [35]. A qualitative discussion of how this 
can be achieved has rcently been provided [36]. Momen- 
tum exchange between the background electron distribu- 
tion and the hole has been made responsible for defor- 
mation of the phase space shape of the hole, with the 
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dynamics of the hole depending sensitively on the shape 
of the background electron distribution. In this Letter we 
present a more quantitative mechanism which is developed 
for electron holes. However, in a similar way it should also 
work for ion holes in the presence of, say, ion conies, which 
have been found in multitude under auroral conditions. 

Mechanism. — Under auroral upward current con- 
ditions the bulk distribution is kind of a non-symmetric 
(downgoing) ring distribution with loss cone (due to the 
presence of the absorbing ionosphere) as shown in Fig- 
ure 1. We assume that some appropriate instability gener- 
ates an electron hole that propagates along the magnetic 
field B at velocity Vh <C Vr much less than the nomi- 
nal electron ring velocity Vr. The question, in which way 
the hole is generated is of secondary importance for the 
purpose of this Letter. The qualitative discussion of [36] 
referred to the two-stream instability [3] as generator of 
the hole. Under the dilute plasma conditions in the up- 
ward current auroral region where uo ce /uo pe ~ 10 — 25, 
the two-stream instability does not grow very fast. Its 
growth rate 7# ~ uob ~ (m e /16 77^)3 ~ 0.03 co pe is 
less than the global electron cyclotron maser growth rate 
7ecm ~ (10- 4 - I0~ 3 )u; ce (cf, e.g. [27-29,40]). However, 
hole generation is not affected by the global maser insta- 
bility which has only a minor effect on the bulk electron 
distribution not causing a substantial energy loss for the 
electrons. Quasilinear flattening of the distribution which 
partially fills in the inner part of the ring (shown as weak 
background in Fig. 1) just where the hole is located in 
phase space is mainly caused by VLF turbulence generated 
under the same conditions [15,35]. However, at fixed v± 
the electron drift velocity Vd might not exceed the electron 
thermal speed v e , in which case the drift cyclotron insta- 
bility (MTSI) takes over. Its growth rate is of the order of 
the lower hybrid frequency j mts ~ {jj lh ~ oo pi which in the 
dilute plasma is also small. Hence, hole formation is a pro- 
cess comparable to or slower than the emission of radiation 
by the global maser instability. In the auroral kilometric 
radiation the latter might thus provide the background 
radiation level while the steep gradients produced in hole 
deformation [36] generate the narrow intense short lived 
drifting emission bands that have been observed. 

The hole-related localised deficiency of electrons on the 
electron background in phase space is centred at the in- 
stantaneous velocity Vh(t) of the hole (Fig. 1). In complete 
analogy to solid state physics it represents a localised pos- 
itive charge Qh in the electron fluid that is attached to the 
hole. The hole moves with velocity Vh(t) on the electron 
fluid. In the presence of an electric field E\\ this charge 
Qh will become accelerated, with its collisionless dynam- 
ics being described by the equation of motion 
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Fig. 1: Schematic of a ring distribution with loss cone. In this 
representation the ring is given a varying phase space density 
with increasing pitch angle ending at a large empty loss cone. 
Also shown is an electron hole at location of parallel velocity 
Vh and for a range of perpendicular velocities v±. The hole 
is assumed to be a straight line initially at constant Vh being 
located between the hot electron ring distribution and the cold 
ion distribution. The latter propagates into opposite direction 
to the electrons at much lower speed. The hole is slow against 
the electrons while its velocity might be comparable to the 
ion speed. Also shown are the components of the nominal ring 
velocity Vr a line of constant v± used in the calculations. 



that holds at every fixed perpendicular velocity v±, as in- 
dicated on the left hand side. The interaction of the hole is 
purely electrostatic, and the electric field is external to the 
hole given by the conditions in the plasma. These are de- 
scribed by the presence of the downward electrons which 
move with respect to the ambient ions a shown in the fixed 
laboratory reference frame in Figure 1. Moreover, the ion 
and electron beams are distributed about homogeneously 
over an area that is large against the extension of the hole 
in any direction. This implies that the electric field can be 
expressed through the current flowing in the plasma as 



dt 



eN 



where v d = (v e \\) - (vq). (2) 



Here N is the plasma density. The drift speed is the differ- 
ence between the average bulk speeds of the electrons and 
ions. In order to find an equation for the parallel velocity 
of the hole, we take the time derivative of Eq. (1), remem- 
bering that the hole charge Qh(vh) itself is a function of 
the hole speed. Defining q = Qh/e and m = Mh/m e this 
yields 



2 Vh _ d\nq(v h ) ( dv h \ 



dt 2 



dv h 



dt J 



q(vh) 



m 



(3) 



Since the hole can move only parallel to the ambient mag- 
netic field, this is a one-dimensional equation of motion 



In deriving this expression we used Eqs. (1-2) in order to 
eliminate the electric field, and we have suppressed the 
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index v± understanding that Eq. (3) holds for every con- 
stant perpendicular velocity. To be able to proceed we 
must determine the charge and mass of the hole. The to- 
tal mass of the hole is given by Mh = m e N t Vh, with N t 
the number density of electrons that are trapped in the 
spatial volume Vh of the hole. This number is assumed to 
be constant during the evolution of the hole, a simplifi- 
cation which neglects any possible exchange of electrons 
between the hole and its environment which determines 
its life time. Hence over the entire evolution of the hole 
m = const does not change. This is not so for the charge 
of the hole Qh = e(N\ v — N t )Vh- Any acceleration of the 
hole in real space corresponds to a displacement in veloc- 
ity space. Hence Qh depends on the fraction in ambient 
number density of ambient electrons of different velocity 
seen by the hole when it moves in velocity space. We need 
to know only the ratio 



q/m = (N\-N t )/N t , 



(4) 



and, hence, the volume Vh of the hole cancels out leaving 
us with the problem to calculate iV|^^ . In order to do 
this we assume, for simplicity, that the background ring 
distribution is a rotated Maxwellian 



F(v\\,v±) 



N 



exp 



r (v-v*n 
i 2*? / • 



(5) 



Here Vr = const is a constant radius in velocity space, 
and v = (v\\ , v_\_) is the velocity vector. (In global electron- 
cyclotron maser theory and simulations as, e.g., in [29], the 
thermal width of the distribution is not important and a 
^-ring distribution is used.) At constant v± the variation 
is in the parallel component (v\\ — Vr\\) 2 in the argument of 
the exponential. The fractional number density of the hole 
is obtained by integration just over the velocity volume of 
the hole. Since this is very small we assume that the hole is 
of rectangular shape of width Avh = \/2cj)h/m e in parallel 
velocity. Then 



N\ Vh = 2irAv h F\ Vh v ± dv ± . 



(6) 



The indication Vh on F means that in F the parallel ve- 
locity v\\ — > Vh is to be replaced by the hole speed. It 
is then easy to show that the coefficient of the second 
term on the left in Eq. (3) is dlng/dv^ = — (vh — )/^e ? 
which is linear in Vh- This fact would enable us to rewrite 
Eq. (3) in terms of the variable x = Vh — Vr\\ . However, 
it is more convenient to introduce dimensionless variables 
r = u pe b,u = v h /v e ,U R = V R \\/v e ,U d = v d /v e ,n = N/N t 
and rewrite Eq. (3) as 

u"-{U R -u){u') 2 = (U d -u) {nCe-^ u «- u ^ - l} , (7) 

where 1 = d/dt, and C for fixed v± is a constant fac- 
tor that is determined from the background distribution 
F. The right-hand side of this equation is linearly pro- 
portional to the average parallel beam velocity U d which 



is the driving force. U d = corresponds to a symmetric 
ring distribution with no loss-cone. The rigth-hand side 
vanishes for U d = u, in which case the equation is solved 
trivially (expressed in decaying error functions in u), and 
it can be shown that an initially finite velocity hole will 
come to rest. Generally one may suspect that asymptoti- 
cally u — > U d . Hence, for a non-symmetric ring distribution 
that, for instance, involves a loss cone, the section of the 
hole corresponding to v± = const is accelerated parallel to 
B until settling near U d . This causes deformation of the 
hole in phase space and is described by the solutions of 
Eq. (7) for the most interesting case U d ^> u. 

The hole is a slowly moving entity of velocity u <^Ur. 
Observation of the displacement of the fine structure in 
the auroral kilometric radiation suggests that holes move 
at velocity Vh < 100 km/s [25] compared to electron beam 
velocities of Vr\\ ~ 10 4 km/s. Hence, in the coefficient of 
the second term on the left in Eq. (7) u can be neglected. 
Moreover, the expression in front of U d — u on the right- 
hand side becomes a constant A = nC exp(— Ur/2). If we 
assume that U d ^> u, which holds for a non-symmetric 
electron distribution like that shown in fig. 1 , the solution 
of Eq. (7) can be found by multiplication with u' and inte- 
grating once with respect to dr. This yields the following 
expression for 



u'(t) = ^/AU d /U R tan \y/AU R U d (T f - r) 



(8) 



with ' 



the final time when the hole comes to rest and 



0. The restriction on the argument of the tangent 



function implies the following restriction on rf. 



Tf <w (2y/AU R U d ) 1 



(9) 



The larger the electron drift velocity U d , the faster the hole 
tends to reach its final velocity. In other words, the faster 
the electron beam moves, the less time it takes for the de- 
formation of the hole in phase space. This is what has been 
expected from the very beginning. Eq. (8) also shows that 
the acceleration of the hole is a monotonically decreasing 
function of time r until the hole arrives at its final state. 
The velocity u(r) of the hole is obtained by integrating 
Eq. (8) with respect to time. In the initial state for times 
r short against rj a solution is found by expanding the 
integrand. Restricting to the first two terms yields 



(10) 



where f = t^/AUrU^. Initially the hole velocity increases 
linearly with time, and the acceleration is proportional to 
\fU dl the root of the electron drift velocity and vanishes for 
U d — > in the case of a symmetric electron ring confirming 
the former conclusion on this case. 

In the more general case when u is of the same order 
as Ur, Eq. (7) cannot be solved analytically. Since this 
case is of lesser interest and would exceed the limits of 
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a Letter, we leave its investigation for a later publication 
(Jaroschek et al., to be submitted). Instead we proceed to 
investigate the variation of the hole velocity with the shape 
of the electron distribution. This is closely related to the 
variation with v± and the drift velocity Ud- To this end we 
study the effect of A = nCexp(— U%/2) — 1 respectively 
C. The constant C is a differential in v± 



n 2Av h 
C = — — — exp 



/7TV e 



(v± ~ Vr±) 2 
2vl 



v±dv^ 



(ii) 



At a phase space section parallel to v\\ all particles in the 
distribution have the same v± — Vr± , and the exponential 
factor in C is unity. We may define the electron drift speed 
and electron temperature at constant v± as the first and 
second moments of the electron distribution function at 
constant v±, respectively. The local drift velocity entering 
into the expressions for the evolution of u{r) is the parallel 
component of the first moment of the electron distribution 
function at constant v±. For a half ring distribution with 
Vm\ > one has 



v±dv± 



(12) 



an expression that decreases with increasing v± . The same 
expression holds for the part of the ring in Fig. 1 that has 
no expression at negative velocities in the loss cone on the 
left of Fig. 1. The decrease in Ud with v± is crucial for the 
evolution of the parallel velocity u{r) of the hole in phase 
space in the presence of a half ring. Since the hole veloc- 
ity and acceleration of the hole in parallel direction are 
proportional to Ud the decrease implies that the largest 
variation in u is for small v±. The velocity space deforma- 
tion decreases with increasing perpendicular velocity v±. 
Hence, the hole attains the largest velocity at small v± 
thus becoming bent in velocity space, as was suggested 
earlier [36] from qualitative considerations. 

If the loss cone is not empty, the effective drift velocity 
Ud = U~£ — U^ the hole experiences is reduced by the 
contribution of the ring at negative velocities 



u d \: 



2%a 



V 2 R 



v±dv± 



(13) 



where a is the fractional density of the half-ring distribu- 
tion with Vr\\ < 0, while the inclusion of the upward ion 
distribution at low v_\_ would increase the effective drift 
speed for the small values of v± to which the ions extend 
in velocity space (see Fig. 1). 

The case Ud < u is not of vital interest. In this case 
u — > 0, and the hole will be about at rest. Such cases may 
refer to standing narrow band emissions in the auroral 
kilometric radiation with the loss-cone about filled. 

Discussion. — The above considerations show that an 
electron phase space hole experiences a particular dynam- 
ics in phase space that depends on the shape of the dis- 
tribution function. This Letter was restricted to the mere 



investigation of the possibility of deformation of an ini- 
tially straight (in velocity space) phase space hole (BGK 
mode) by momentum exchange with the ambient electron 
distribution. We argued that momentum exchange causes 
differential acceleration of the hole different for different 
v±. 

The evolution of phase space holes has so far been con- 
sidered in theory and simulations only for Maxwellian dis- 
tributions in which case their phase space dynamics is 
simple (see, e.g., [23,33]). Observations in space, for in- 
stance under conditions in the auroral kilometric radiation 
source and also in collisionless shocks, suggest that phase 
space holes evolve when the phase space particle distribu- 
tions deviate strongly from Maxwellian shape [6,10-12]. In 
particular, the combination of an electrostatic field along 
the magnetic field and a magnetic mirror geometry trans- 
form a beam distribution into a ring distribution [5] with 
loss cone (colloquially called a horseshoe). Such distribu- 
tions are the rule in the aurora and are also expected in 
super-critical quasi-perpendicular collisionless shocks [37] . 
In their presence a hole should undergo deformation in 
phase space of the kind described in this Letter. Holes have 
also been predicted in relation to reconnection in collision- 
less current sheets [7] where their signatures might have 
been observed in situ in space [34,39]. 

The phase space deformation might not be of over- 
whelming importance for the dynamics of the plasma even 
though its importance for dissipation has not yet been 
investigated. However, as suggested in [36] it should be 
of crucial importance for involving electron holes into 
the emission of electromagnetic radiation by the electron- 
cyclotron maser mechanism acting in the auroral regions 
of planetary magnetospheres and in a variety of astrophys- 
ical objects (cf., e.g., [1]). This requires the generation of 
sharp positive phase space gradients dF(v\\\,v±)/dv± > 
on the electron distribution function in the perpendicular 
velocity component. Phase space holes do not originally 
possess such gradients and therefore should not become 
involved into generation of radiation other than contribut- 
ing to electron heating and acceleration of electron beams 
which on their own might secondarily generate plasma 
waves away from the hole and produce second harmonic 
plasma radiation with frequency 2co pe . This mechanism is 
viable, for instance, at collisionless shocks and might be 
responsible for the so-called backbone radiation in type II 
radio bursts observed in the solar corona and interplane- 
tary space. 

However, for electron holes to become directly involved 
into the electron cyclotron maser mechanism, generation 
of perpendicular gradients is absolutely necessary. This is 
provided by the phase space deformation mechanism pro- 
posed in this Letter, which bends the hole in phase space 
and transforms its natural parallel phase space gradient 
dF(v\\,v±)/dv\\ into a perpendicular gradient. Since 
dF(v\\,v±)/dv\\ 7^ is steep and increases when the hole 
is accelerated into the bulk of the distribution [36], the 
new dF(v\\\,v±)/dv± is also very steep and will readily 
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contribute to maser emission. Its frequency maps the lo- 
cal electron cyclotron frequency or its lower harmonics. 
We are not going here into the details of the radiation 
process as this has been described in the literature for 
the global electron distribution (see e.g. [18,19,27]). Since 
the hole possesses two boundaries, one of them causes a 
positive, the other a negative phase space gradient. It has 
been argued [36] that, theoretically, these gradients corre- 
spond to emission and absorption separated by the whole 
width. Current instrumental resolution and possibly nat- 
ural line broadenings do not allow for a discrimination of 
these emissions and absorptions, however. Whistler saucer 
emissions from holes might be another indication of evo- 
lution of positive perpendicular phase space gradients. 

The emission frequency changes when the hole is dis- 
placed in space along the magnetic field. Because of the 
steepness of the gradient, the emission is very narrow 
band. From its spectral displacement the electron hole 
speed can be determined while from its spectral width 
properties of the hole can be inferred [36]. In the auroral 
region the hole speed is found to be of the order of < 100 
km/s, which is comparable to the ion- acoustic speed be- 
ing sufficiently far below the electron thermal velocity to 
justify the model of a slow hole. However, the measured 
spectral speeds reflect the bulk velocity of the hole. This is 
determined by integrating over the hole distribution along 
the bent phase space shape of the hole. It can be much 
less than the differential speeds of its fastest moving small 
v± phase space sections. 

We also note that ion holes behave in a similar way ex- 
periencing bending in phase space if the ion distribution 
differs from a Maxwellian. This is the case in regions where 
ion conies are generated in the presence of field aligned 
electric potentials and mirror geometries for upward go- 
ing ions. The evolving steep perpendicular velocity space 
gradient may excite ion cyclotron waves from ion holes. 
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